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The top quark, discovered in 1995 by both Tevatron 
experiments [l[ , is the only known fermion with a mass of 
the order of the electroweak breaking scale. This suggests 
that it may play a special role in new physics. A detailed 
investigation of the production mechanism of top quarks 
will give insights into whether top quarks are produced 
via new physics processes. 

In this Letter we present two analyses studying the 
forward-backward charge asymmetry of top quark pairs 
in pp collisions at y/s = 1.96 TeV at the Fermilab Teva- 
tron. In standard model QCD a charge asymmetry Ac 
arises in next-to-leading order ttX production. Because 
the strong interaction is invariant under charge conjuga- 
tion Aq is equivalent to a forward-backward asymmetry 
Afb ■ Recent calculations predict a slightly positive total 
Afb = 5.0 ± 1.5% in the Tevatron pp rest frame @, 
with the theoretical uncertainty driven by the size of cor- 
rections at higher orders. This small total Afb combines 
a positive asymmetry from the interference of the Born 
and virtual (box) corrections (ti) with a negative asym- 
metry from interference of initial and final state radiation 
amplitudes (ti + g) Q . 

While the total A-pB value expected by the standard 
model is hardly measurable at the presently achievable 
precision, we are sensitive to large Ayb values (of order 
±30%) predicted in some models with new physics, e.g. 
Z'-like states with parity violating couplings [5j and the- 
ories with chiral color [2J, |6j . In contrast to searches for 
heavyresonances in the spectrum of the mass of the top 
pair [7[, a measurement of Apb is sensitive to both nar- 
row and broad resonances. In addition, the presence of 
a massive gluon may be visible in the asymmetry even 
above the collision energy due to interference with the 
standard model gluon. 

Since a longitudinal boost changes the top quark di- 
rection, Afb is frame dependent. Undetected collincar 
gluon radiation makes the fundamental initial parton 
frame experimentally inacessible. However, the tt and 
the frame are experimentally accessible and according 
to [2| the Afb values in the pp frame are predicted to be 
reduced by w 30% relative to the tt frame. 

We present here the first measurement of the top quark 
production A-pBi fully corrected to the parton level, in 
both the pp and tt frames. Correction to the intrinsic 
parton value allows direct comparison to theoretical pre- 
diction, and measurements in two frames probe the con- 
sistency and the frame dependence of the effect. A recent 
study [8[ measures a quantity which is related to the ti 
frame asymmetry but is uncorrected for acceptance and 
resolution effects. The result (12±8±1)% is larger than 
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expected, within errors, but difficult to interpret. 

We use 1.9 fb _1 of pp collision data recorded by the 
CDF II detector. The detector is a forward-backward 
symmetric system consisting of a magnetic spectrometer 
surrounded by projective calorimeters and muon detec- 
tors 0. Charged track reconstruction in a 1.4 T axial 
field uses a large open cell drift chamber and silicon mi- 
crostrip detectors for displaced secondary vertex detec- 
tion. We use coordinates where <fi is the azimuthal angle 
and 6 is the polar angle with respect to the proton beam 
z-axis. Transverse energy is Et = E sin 6, the rapidity is 
Y = i \n[(E + p z )/(E ~ p z )], and the pseudorapidity is 
r] = -ln[tan(0/2)]. 

We collect a sample of candidate events in the lep- 
ton+jets topology tt (W + b)(W~b) -> (qq'b)(£- D e b) 
[l(|, where one W-boson decays leptonically and the 
other hadronically, by triggering on a central (\rj\ < 1.0) 
electron with Et >18 GeV or central muon with trans- 
verse momentum px >18 GeV/c. After offline recon- 
struction we select events with an isolated electron with 
Et _ 20 GeV or muon with pt > 20 GeV/c, missing 
transverse energy &t _ 20 GeV |ll[ consistent with a 
neutrino from W decay, and at least four hadronic jets 
with |?7| < 2.0 and Et > 20 GeV. Jets are clustered 
in fixed cones of radius AR = ^(Ar/) 2 + (A0) 2 < 0.4 
and jet energies are corrected to parton-level values [13] • 
At least one jet must be ^-tagged, i.e., contain a recon- 
structed secondary vertex consistent with the decay of 
a bottom hadron in the jet [Hj]. We find 484 candidate 
events. 

The expected ti signal is studied using the pythia, 
herwig, and MC@NLO event generators [lj] and a full 
detector simulation [l5T ] . The top quark mass is set equal 
to M t = 175 GeV/c 2 . The rates and kinematics of back- 
ground processes are well modeled with simulation and 
data control samples [l6| which will be discussed later. 
We expect a total of 87±23 background events, leaving a 
ti signal of 397 ± 32 events, consistent with our previous 
cross-section measurement of 8.2 ± 1.0 pb [17] . 

Mass constraints on the W bosons from top quark de- 
cay fix the jet parton assignment and allow complete 
reconstruction of the ti kinematics. For the pp frame 
analysis we use the algorithm employed in the top quark 
mass measurement of Ref. JLj|. Measured jet energies 
float within expected resolutions, 6-tagged jets are taken 
as fragmented 6-quarks, both W boson masses M(tv) 
and M(qq') are constrained to 80.4 GeV/c 2 , and the 
top quark mass is constrained to 175 GeV/c 2 . For 
the ti frame analysis we use the technique described in 
Refs. [lj| which employs constraints on the W boson 
masses, the reconstructed t-i mass difference (but not 
M t ), the total transverse energy, and the ^-likelihood 
of the jets [20j ]. In simulated ti samples the two pro- 
cedures resolve the top direction with similar accuracy. 
The resolution on the direction of the hadronically de- 
caying top quark th, expressed in terms of rapidity, is 
ay(th) — 0.29. The leptonically decaying top quark 
system ti, which includes the indirectly measured neu- 
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TABLE I: Measured asymmetries in large simulated tt samples. 



Generator 


Parton-level Ap 1 ^ 


Reconstructed A 7 ^ 


Parton-level Ap B 


Reconstructed Ap B 


Pythia 


0.000 ±0.003 


-0.007 ±0.006 


0.000 ± 0.001 


-0.005 ± 0.003 


Herwig 


0.000 ±0.006 


-0.013 ±0.012 


-0.003 ± 0.002 


-0.003 ± 0.006 


MC@NLO 


0.038 ± 0.002 


0.015 ±0.016 


0.049 ± 0.002 


0.017 ±0.007 



trino, has ayite) — 0.46 and significant non-Gaussian 
tails (15%). 

We measure the direction of the top quark in the pp 
center-of-mass frame using the cosine of the polar angle 
between the hadronic top quark and the proton beam, 
cosckp. The sign of the t^ electric charge is opposite 
that of Qe, the leptonic charge observed in the i; de- 
cay. Assuming CP invariance, we find one top quark 
angle cos 9 — —Qe • cos a p in each event and calculate the 
asymmetry in the pp center of mass frame plj 

a pp _ ^(cos 9 > 0) - N(cos6 < 0) 
FB ~ iV(cos6> > 0) + N(cos9 < 0)' ( ' 

This technique has the simplicity of relying only on the 
hadronic top quark reconstruction, but has the drawback 
of measuring asymmetries which are diluted by 30% com- 
pared to the tt frame. 

The tt rest frame measurement exploits the Lorentz 
invariant difference between the t and i rapidities Y t and 
Yf . We use the reconstructed rapidity of tt and th in 
each event, assume CP invariance, and determine AY — 
Y t — Y[ = Qe ■ (Y t( — Y th ) from which we calculate the 
asymmetry in the approximate (LO) tt rest frame (22j 

A tt _ N(AY > 0) - N(AY < 0) 
FB ~ N(AY > 0) + N(AY < 0) ' () 

To connect this with other asymmetry measurements, we 
note that in the case of ideal resolution Ap B reproduces 
the asymmetry measured in the equivalent Collins-Soper 
frame 23]. While it is sensitive to the larger tt frame 
asymmetry, Ay combines the uncertainties of both quark 
reconstructions, including the neutrino-related complica- 
tions of the tt quark system. 

The expected measurement performance is evaluated 
using simulated samples. In Table[I]we compare asymme- 
tries found after selection and reconstruction to parton- 
level asymetries calculated using perfect acceptance and 
resolution. The uncertainties reflect the simulation 
statistics. With the parton-shower generators pythia 
and herwig we see no intrinsic charge asymmetry at the 
parton-level, as expected, and verify that any forward- 
backward bias from selection and reconstruction is small. 
With the MC@nlo generator, which includes the small 
QCD-induced charge asymmetry, we find parton-level 
values consistent with theoretical expectation in mag- 
nitude and the level of frame dependence. With large 
statistics the measured values are sensitive to the small 
asymmetry, but diluted by acceptance and reconstruc- 
tion effects. Dilution corrections, as well as the expected 



sensitivity in our finite dataset, are discussed later. The 
calibration of the simulation to the physical detector ge- 
ometry and acceptance has been checked in studies of 
electroweak processes [24| , For example, the leptonic 
charge asymmetry in — * V^v agrees with our sim- 
ulated physics and detector model within the statistical 
uncertainty of w 0.004. 

The Afb measured in data must be corrected for back- 
ground contributions which include asymmetric weak 
processes. W + jets events with tagged heavy flavor 
(W+hf) or mistagged light partons (W+lf) are modeled 
using ALPGEN [25| interfaced to pythia parton shower- 
ing, along with ^-tagging and mistagging rates parame- 
terized from jet data. Small electroweak backgrounds 
(EW), WW, WZ, ZZ and single-top, are modeled with 
pythia and with MADEVENT |26| respectively. The non- 
14 7 (QCD) electron background is studied using data 
events with five jets where one jet models a misrecon- 
structed electron; the same sample is used for non-14 7 
muons after re-weighting the lepton acceptance. The 
background levels and asymmetries expected in the two 
analyses are shown in Table [ill The combined results 
are listed in the last row. 



TABLE II: Backgrounds, estimated number of events N e: 



and their effective asymmetries. 


Process 

W+hf 

W+lf 

EW 

QCD 


■^V exp 

37 ± 10 
20 ±5 
12 ± 1 
18 ±16 


APP 

-0.087 ±0.005 
0.044 ± 0.008 
-0.022 ±0.014 
-0.008 ± 0.004 


Att 

— 0.045 ± 0.003 
-0.006 ±0.015 
-0.015 ±0.044 
0.006 ±0.010 


Total 


87 ±23 


-0.053 ± 0.004 


-0.021 ± 0.007 



Fig. [T] shows the measured distributions of cos 9 and 
AY in the 484 6-tagged tt candidates, along with pre- 
dictions based on simulated ti events from the mc@nlo 
generator in combination with our non-tt background 
models. The measured asymmetries are displayed in Ta- 
ble IIIII The background-corrected values, derived by 
subtracting the composite model shape bin-by-bin, show 
a positive asymmetry which is larger than but consis- 
tent with the MC@NLO predictions within uncertainties. 
Our background-corrected Ap B , although measured in a 
slightly different visible phase space, is very consistent 
with the measurement from the D0 collaboration 
Subdividing the data by lepton types and lepton charges 
shows a consistent positive asymmetry across all samples. 

To study the two contributions tt and tt + g with dif- 
ferent expected sign in A^b, we split our data sample 
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FIG. 1: The two top quark production angle variables, cos# for the pp frame and AY for the 
tt frame. The solid line is the prediction for tt with mcQnlo model of the QCD induced charge 
asymmetry and er t j = 8.2 pb, plus the expected non-tt backgrounds. The dashed curve shows the 
prediction when tt is reweighted according to the form 1 + Afb cos a using measured values of Afb- 



TABLE III: Predicted (mcQnlo + non-t£) and measured asymmetries without further corrections. 





Nj e ts 


Predicted A% 


Measured Apg 


Predicted A" B 


Measured A P * B 


Recontructed Data 


> 4 


0.001 ±0.010 


0.099 ± 0.045 


0.010 ±0.007 


0.087 ± 0.045 


Background Subtracted 


> 4 


0.015 ±0.016 


0.130 ±0.055 


0.017 ±0.007 


0.119 ±0.064 




= 4 


0.032 ±0.018 


0.120 ±0.064 


0.038 ± 0.008 


0.132 ±0.075 




> 5 


-0.027 ±0.032 


0.160 ±0.109 


-0.033 ±0.012 


0.079 ±0.123 



into events without any additional hard jet (Nj e t s = 4, 
85% it) and events with at least one additional hard jet 
(Nj e ts > 5, 53% tt + g). Our background-corrected A^ B 
and Ap B values for this study are presented in Table IHIl 
The Nj e t s dependence is not as strong as seen in [H, but 
the limited statistics does not allow a firm conclusion. 

The distributions in Fig. [1] are distorted from their 
true parton-level shapes by acceptance bias and recon- 
struction errors. We use a matrix inversion technique to 
derive the parton-level distributions and tt asymmetries. 
If an event in bin j at parton-level is collected with effi- 
ciency Ej and migrates to bin i at the measurement level 
with probability Sy, the bin-by-bin parton-level distri- 
butions Pj can be found from the background-corrected 
data distributions Di by the inverse transformation 

Pj = s Tl • V • n - ( 3 ) 

We simplify each distribution to four bins, with two bins 
on either side of the cross-over at cos 9 = AY = 0. The 
efficiencies and migration matrix are derived by com- 
paring the parton and reconstructed level quantities us- 
ing the zero asymmetry pythia tt simulations. In the 
cos 9 (/AY) analysis roughly 13% (25%) of events change 
signs, but the matrix is symmetric within uncertainties. 
The symmetry of the matrix, which follows from the 
forward-backward symmetry of the detector, ensures that 



the inversion is insensitive to small errors in the modeling 
of the migration parameters. 

The expected performance of the complete calculation 
is evaluated with simulated samples. Sensitivity to the 
asymmetry model is studied using pythia samples that 
have been reweighted in the top quark production angle 
for a range of possible asymmetry functions and magni- 
tudes varying between 0.0 and 0.30. Sensitivity to the 
QCD-induccd asymmetry is studied with MC@nlo. The 
effect of extra jets is studied with exclusive it + par- 
ton and it + 1 parton samples made with the ALPGEN 
generator. Each sample was reconstructed, measured, 
and propagated back to the parton level with the proce- 
dures described above. For all conditions the procedure 
returns mean values within 0.02 of the true value. The 
predicted statistical precisions in our 1.9 fb _1 dataset are 
5A% = 0.09 and 5A% = 0.13. 

Additional sources of uncertainty are evaluated using 
simulated samples with reasonable variations on the as- 
sumptions for background shape and normalization, sig- 
nal shapes, the top quark mass, the parton distribution 
functions, the amount of initial and final state gluon ra- 
diation, and the calorimeter energy scale. The largest 
uncertainty in the A^ P B analysis is the background nor- 
malization and the largest in the A^ B analysis is the 
AF shape modeling, being roughly SAfb — 0.02 in each. 
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The total systematic uncertainty is 5A^b = 0.04 for both 
techniques. 

Applying our algorithm to the inclusive background- 
subtracted distributions in Fig. [TJ we find parton level 
asymmetries of A% = 0.17 ± 0.07 ± 0.04 and A% = 
0.24 ± 0.13 ± 0.04, where the uncertainties are statistical 
and systematic, respectively. In Fig. [TJ the dashed lines 
show that the data are in good agreement with models 
derived by reweighting the generated top quark produc- 
tion angle a in the symmetric pythia sample with form 
1 + A FB cos a using the measured A FB . 

In conclusion, we have measured a forward-backward 
and (equivalent) charge asymmetry in a strong process at 
high energy using reconstructed tt events produced in pp 
collisions at y/s = 1.96 TeV. We find forward-backward 
parton-level asymmetries of A F P B = 0.17 ± 0.08 in the pp 
frame and A FB — 0.24 ±0.14 in the ti frame. Our results 
show the expected frame dependence, are consistent (< 
2a) with the small (~ 0.05) charge asymmetry expected 
from QCD, and they disfavor exotic sources of top quark 
production with significant negative Ayb values [2| . 
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